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, Abstract. Density profiles of cosmological virialized systems, or dark 

CN ' halos, have recently attracted much attention. I first present a brief 

historical review of numerical simulations to quantify the halo density 
profiles. Then I describe the latest results on the universal density profile 
and their observational confrontation. Finally I discuss a clustering model 
of those halos with particular emphasis on the cosmological light-cone 
effect. 
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' 1. Introduction 

o . 

■ The key assumption underlying the standard picture of structure formation is 

that the luminous objects form in a gravitational potential of dark matter ha- 



X 



los. Therefore, a detailed description of halo density profiles as well as of their 
clustering properties is the most basic step toward constructing the formation 



O ' and evolution of galaxies and clusters. 



^ ' of dark halos is two- fold: 



More specifically, the importance of the detailed studies of density profiles 



(i) theoretical interest; 

What is the final (quasi-)equilibrium state of cosmological self-gravitating 



' systems (as long as the energy dissipation is neglected) ? One may easily 

think of two quite distinct, but equally plausible, possibilities; 

(A) the systems reach a certain universal distribution which is indepen- 
dent of the cosmological initial condition. 

(B) the systems somehow keep the memory of the cosmological initial 
condition even at the highly nonlinear regime. 

The singular isothermal sphere: 

may be a reasonable possibility along the line of the idea (A). As a matter 
of fact, quite often I meet people who argue that the idea (B) is unlikely 
because of the strong nonlinear nature of the gravitation. In such an 
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occasion, I present an example of the well-known stable clustering solution 
for the nonlinear two-point correlation function (Davis &: Peebles 1977): 

Prasssik) oc A;" ^ ema.s(r)cxr-3("+3)/("+5). (2) 

This provides a good specific case that the cosmological initial condition 
is not erased and imprinted even in the strongly nonlinear behavior as is 
well confirmed by later numerical simulations (e.g., Suto 1993; Suginohara 
et al. 2001). Of course the answer to the final state of cosmological self- 
gravitating systems may not be unique since it should really depend on 
the specifics and scales of the systems under consideration. For instance, 
it is clearly hopeless to extract any meaningful cosmological information 
from the precise data on the orbit of the earth around the Sun; all the 
initial memory should have been lost due to the strongly nonlinear and 
chaotic nature of the gravitation. Nevertheless the final state of dark 
matter halos corresponding to galaxy- and cluster-scales is a well-defined 
problem which may be reliably answered with the current high-resolution 
numerical simulations independently of the physical intuition. 

(ii) practical application; 

Whether or not the density profiles of dark halos keep the cosmological 
initial memory, the quantitative (empirical) prediction of the profile for 
a given set of cosmological parameters has several profound astrophysical 
implications including the rotation curve of spiral galaxies, reconstruction 
of the mass distribution from the weak-lensing, and X-ray and SZ observa- 
tions of galaxy clusters. In particular, the confrontation of those testable 
predictions against the accurate observational data may even challenge the 
cold dark matter paradigm itself. 

In this article, first I will review the summary of the past studies of the 
density profiles of dark matter halos, and then present some applications of 
those results. 

2. Previous work on the density profiles of dark matter halos 

The study of the density profiles of cosmological self-gravitating systems or dark 
halos has a long history. Table 2 is my personal (and thus incomplete) summary 
of some past work that I recognize (or at least I have really read those papers). 

Hoffman & Shaham (1985) argued, on the basis of the secondary infall 
model (Gunn & Gott 1972), that the density profile around peaks is given by 
p oc 7--3(n.-l-3)/(n-i-4) -^yj^gi-g ^ ^ d\n P{k)/d\nk is the spectral index of the density 
fluctuation field. The resulting power index — 3(n -|- 3)/(n -|- 4) is close to, but 
slightly differ from, that for the stable clustering solution for the nonlinear two- 
point correlation function — 3(n-|-3)/(n-|-5). In any case this result implies that 
the density profile of dark halos does keep the initial memory. 

Many subsequent numerical simulations seem to have confirmed the pre- 
diction of Hoffman &; Shaham (1985). While some authors have reported the 
existence of a kind of universal profile, they were not able to specify its func- 
tional form. In this respect, the proposal of the specific universal density profile 
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by Navarro, Frenk & White (1995, 1996, 1997) was really a breakthrough. They 
suggested that all simulated density profiles can be well fitted to the following 
simple model (now referred to as the NFW profile): 

(r/rs)(l + r/rs)^ 

by an appropriate choice of the scaling radius Vs = rs{M) as a function of the 
halo mass M. Considering the previous conclusion from numerical simulations, 
their claim was quite surprising and, in some sense, contradictory. It turned 
out that most simulations prior to their discovery did not have sufficient mass- 
resolution to reliably represent the deep central gravitational potential, which 
gives rise to an artificial central core. This presents a good example that the 
quantitative difference leads to the qualitatively different interpretation of the 
physical phenomena; now the density profile of dark halos seems to lose the 
initial memory, and thus the answer to the final (quasi-) equilibrium state of 
cosmological self-gravitating systems changes even in a qualitative way; from 
the possibility (A) to (B) mentioned in the previous section ! 



Table 1. The stable solution and the density profile around peaks. 

n Stable OC r-3("+3)/{n+5) p ^ ^-3(n+3)/(n+4) 

~o :p''^ 

-2 r'^ r~^-^ 



If this change of the qualitative conclusion (from Hoffman Sz Shaham to 

NFW) is really ascribed to the mass-resolution of the numerical simulations, 
one should naturally wonder whether or not current numerical simulations are 
sufficiently reliable for the present problem. This motivated many people to 
examine the convergence of the profile, in particular, its inner power-law index, 
with significantly higher-resolution simulations. Fukushige &; Makino (1997) 
are the first to claim that the inner slope of density halos is much steeper than 
the NFW value. Their claim was confirmed by a series of systematic studies by 
Moore et al. (1998, 1999), and the current consensus among most (but not all) 
numerical simulators in this field is the converged profile is given by 

p{i^) ^ 7 — ; — r-r. ; — — (4) 

(r/rs)°(l + r/rs)3-° ^ ^ 

with a ~ 1.5 rather than the NFW value, a = 1, (e.g., Fukushige &: Makino 

2001) . In most cases this difference is fairly minor and the NFW profile still 
provides a good empirical approximation to the dark matter halo (Jing & Suto 

2002) . 

I just reproduce two figures from my collaborative work (Jing &: Suto 2000) 
just to illustrate the validity of the universal density profile (eq.[4]) despite the 
fact that their projected images look very diverse rather than universal. It is 
hard to imagine the regularity of the averaged profiles plotted in Figure 2 behind 
the apparent diversity quite visible in Figure 1. 
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Table 2. Past work on the density profiles of dark matter halos. 



1970 Peebles 

N = 300 simulation to reproduce the profile of the Coma cluster 
1977 Gunn 

predicted p oc r~^/^ on the basis of secondary infall model 

1985 Hoffman &: Shaham 

density profiles around the peak: p oc y.-3{"+3)/(n+4) 

1986 Quinn, Salmon Sz Zurek 

confirmed the Hoffman & Shaham prediction (N ~ 1000 simulation) 

1987 West, Dekel & Oemler 

suggested of the existence of a kind of universal density profile 
with N ~ 4000 simulations. 

1988 Frenk, White, Davis k Efstathiou 
reproduced the galactic rotation curve 

in the SCDM model with N = 32^ simulations 
1990 Hernquist 

proposed an analytic model for density profiles of elliptical galaxies 

p{r) = (Ma/27r)r-^(r + a)~^ 

1994 Crone, Evrard, & Richstone 

suggested a central cusp (not a core) with A'^ = 64^ simulations 

1995 Navarro, Frenk & White 

discovered a universality in the halo density profile 
from SPH simulations in SCDM, A'^ ~ 6000 per halo 

pir) = 7500p(r/0.2r2oo)-Hl + '^/0-2r2oo)"^ 

1996 Navarro, Frenk & White 

19 halos with N{< ryir) = 5000 - 10000 in SCDM 
suggested a specific form for the universal density profile 

1997 Fukushige k Makino 

found a steeper inner profile with N = 786400 simulations 

1997 Navarro, Frenk & White 

universal density profile in a variety of cosmological models 

1998 Syer & White 

repeated merger + tidal disruption model — p oc 7--3{n+3)/(n+5) 

1998 Moore et al. 

found a steeper inner profile p oc r~^-^ 

consistent with the previous finding of Fukushige k Makino 

1999 Moore et al. 

suggested that p{r) oc (r/rs)~^-^[l + {r/rs)^'^]~^ 

2000 Jing k Suto 

dependence of the inner slope on the halo mass (controvertial) 
2002 Jing k Suto 

triaxial modeling of the universal dark matter density profiles 
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Figure 1. Snapshots of the simulated halos at z = 0. Top, center 
and bottom panels display the halos of galaxy, group and cluster masses 
(for four different realizations from left to right), respectively. The size 
of each panel corresponds to 2rvir of each halo (Jing & Suto 2000). 




0.01 0.1 



0.01 0.1 



0.01 0.1 



Figure 2. Spherically- averaged radial density profiles of the simu- 
lated halos of galaxy (left), group {middle), and cluster (right) masses. 
The solid and dotted curves represent fits of a = 1.5 and a = 1 respec- 
tively. For reference, we also show p[r) cx r~^ and r~^'^ in dashed and 
sohd lines (Jing & Suto 2000). 
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3. Observational confrontation of the universal density profile 
3.1. Mass profile and rotation curve 

The advantage of the NFW profile is that it is completely specified for a given 
of cosmological parameters. Therefore this leads to a variety of testable quanti- 
tative predictions. In fact, one can write down all the necessary formulae in the 
following concise manner: 



ScPcrit 



(r/rs)(l + r/rs)2' 
where Pcrit(z) is the critical density of the universe at z 



(5) 



PcM = + « 1.8 X 10^1(1 + zfh'M^/Mpc\ (6) 

The characteristic halo density contrast 6c, and the scaling radius are written 
in terms of the concentration parameter as 



Ayir^^O 



3 ln(l + Cvir) - Cvir/ (1 + Cyir) ' 



rJM) . 



Cvir 



(7) 

(8) 



Finally the virial radius ryir = ryir{M, z) and the concentration parameter are 
given as 



3M 



47rAvirf^0Pcrit 



1/3 



1.69 /Avir^^O 



l + z\ 187r2 



M 



1/3 



Cvir(M, Z) 



187r2j^(z)-0-7 (Ao = 0) 
18TT^n{z)-^-^ (Ao = 1 - ^0 
M 



1 + z Vl.5 X 1013/i-iMp 



-0.13 



(9) 
(10) 

(11) 



The expression for Cvir is the fitting formula in LCDM (Lambda Cold Dark 
Matter) by Bullock et al. (2001). 

Then the corresponding mass profile and the circular velocity are computed 
in a straightforward manner as 



M(r) = 47r 
and 

VJr) 



^cPcTitr"^ 



(r/rs)(l + r/rs)- 



;dr = 4:TrScPcritr, 



In 1 + 



r + rs 



'GM{r) 



1 47r^cPcritrf 



- In 1 + - 



r + Ts 



(12) 



(13) 



The above results explicitly show that one can in principle compare the universal 
density profile with the observational data without any unknown free parameter. 
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This is why the reported disagreement with the mass profile reconstructed from 
weak/strong lensing (Tyson et al. 1998) and the rotation curves of dwarf/low 
surface brightness galaxies (e.g., Moore et al. 1999; de Blok et al. 2001), if real, 
cannot be easily accommodated in the context of the standard paradigm of cold 
dark matter. 

3.2. X-ray gas density profile from the universal density profile of 
dark halos: origin of the isothermal /? model ? 

Another important application of the universal density profile of dark halos may 
be found in X-ray emission from galaxy clusters. A simple but realistic approx- 
imation is the isothermal gas distribution in hydrostatic equilibrium embedded 
in the NFW halo gravitational potential. Actually the hydrostatic equilibrium 
equation: 



Pgas, (14) 



dr 

combined with the NFW profile and equation of state 

Pgas ~ '^gas^^gas ~ ^^gas (1^) 

can be analytically integrated (Makino, Sasaki &; Suto 1998; Suto, Sasaki k, 
Makino 1998) as 

Pgas(r) = Pgo (l + ^) ""^^'^ , B = ^47r<5ePcritr^ (16) 

\ Ts J rK,J- gas 

More surprising is the fact that the above model is very close to the conventional 
/? model. Indeed we found the following fit as illustrated in Figure 3: 

^""^'^ [r+FS^)¥^' ^^^^ 

A « -0.1786 + 0.982 = -0.013B + 0.982 (6 = 25/27), 
rceff ~ 0.22rs, /?eff ~ 0.96 = 0.067S. 



As far as I understand, this work is the first analytical model for the cluster 
gas distribution obtained by consistently solving the Poisson equation for the 
underlying halo mass distribution. The fact that the resulting gas profile is 
quite similar to the isothermal /3 model, widely used as a good empirical model, 
is encouraging. On the other hand, the expected core sizes of clusters seem to 
be significantly smaller than the observational ones (Fig. 4). In fact, rotation 
curves of dwarf/LSB galaxies (Moore et al. 1999; de Blok et al. 2001) and the 
mass profile reconstructed for the galaxy chistcr CL0024 (Tyson et al. 1998) 
also indicate the presence of the core, rather than cusp, in the halo density 
profile. While this is still controversial at this point, Spergel h Steinhard (2000) 
among others even proposed to abandon the collisionless nature of dark matter 
to reconcile the reported discrepancy. 
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Figure 3. Gas density profile (solid 
lines) expected from the universal den- 
sity profile of dark matter halo (dashed 
line) for b = 1/3, 1, and 5/3. For 
comparison, the best-fit /3-models with 
/3 = 0.96 and Vc = 0.22rs are plotted in 
dotted lines (Makino, Sasaki, & Suto 
1998). 



Figure 4. Predicted properties of gas 
density distribution as functions of the 
X-ray cluster gas temperature Tx = 
Tyir{rs)- (a) Ratio of ryir and rc(~ 
0.22rs); (b) ratio of the virial temper- 
atures at and at r^^^j.; (c) predicted 
sizes of the effective core radius com- 
pared with the observed cluster data. 
Solid, dotted and dashed lines indicate 
SCDM, OCDM, and LCDM models, 
respectively. 



Tx - T,„(r,) [keV] 



3.3. CoUisional dark matter 

Since the cold dark matter model is known to successfully reproduce the large- 
scale structure of the universe (r ^ IMpc), the possible collision of dark matter 
should be effective only at scales below IMpc. The required collision cross 
section a for a dark matter particle of mass m may be obtained by assuming 

(mn) ( —] I ^ 1 for mn ^ pc duster and I IMpc. (18) 

This may be rewritten as 

. 1.6/.--cmVg ( { (19) 

Km J \Pc,cluster ) \ < / 

or equivalently as 

. » 2.5 X 10-«/.-'an^ {-^) { '^^] ( 'J^] . (20) 
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The corresponding mean collision time-scale is 

I q /10()(Jkm/s\ / / \ , 1 
At « - « loM — h-^year. (21) 

The above rough estimate is encouraging because in the Hubble time collisional 
dark matter with the cross section of an order a/m ^ Icm^/g will affect the 
central part of dark halos only while keeping the larger-scale behavior intact. 

More detailed study of the effect of collision should be made again with 
numerical simulations. Yoshida et al. (1999) and Dave et al. (2001) conducted 
such simulations and found that dark matter models with a/m ^ Icui^/g lead 
to the formation of the softened central core, rather than cusp, in cluster-sized 
halos. The models, however, simultaneously produce too spherical halos to be 
compatible with the gravitational lensing observations (Miralda-Escude 2002). 
In addition, there are other indications to favor the presence of halos on the 
basis of the gravitational lensing statistics (Molikawa & Hattori 2001; Oguri et 
al. 2001; Takahashi & Chiba 2001; Li & Ostriker 2001). Furthermore, CL0024 
seems to be a very complicated system and thus the interpretation of its density 
profile is not straightforward (Czoske et al. 2002). 



4. Beyond spherical description of dark halos 

Actually it is rather surprising that the fairly accurate scaling relation applies 
after the spherical average despite the fact that the departure from the spher- 
ical symmetry is quite visible in almost all simulated halos (Fig. 1). A more 
realistic modeling of dark matter halos beyond the spherical approximation is 
important in understanding various observed properties of galaxy clusters and 
non-linear clustering (especially the high-order clustering statistics) of dark mat- 
ter in general. In particular, the non-sphericity of dark halos is supposed to play 
a central role in the X-ray morphologies of clusters, in the cosmological parame- 
ter determination via the Sunyaev-Zel'dovich effect and in the prediction of the 
cluster weak lensing and the gravitational arc statistics (Bartelmann et al. 1998; 
Meneghetti et al. 2000, 2001; Molikawa & Hattori 2001, Oguri et al 2001). 

Recently Jing &: Suto (2002) presented a detailed non-spherical modeling 
of dark matter halos on the basis of a combined analysis of the high-resolution 
halo simulations (12 halos with N ~ 10^ particles within their virial radius) 
and the large cosmological simulations (5 realizations with N = 512^ particles 
in a 100/i~^Mpc boxsize). The density profiles of those simulated halos are well 
approximated by a sequence of the concentric triaxial distribution with their 
axis directions being fairly aligned. They characterize the triaxial model quan- 
titatively by generalizing the universal density profile which has previously been 
discussed only in the framework of the spherical model, and obtain a scries of 
practically useful fitting formulae in applying the triaxial model; the mass and 
redshift dependence of the axis ratio, the mean of the concentration parameter, 
and the probability distribution functions of the the axis ratio and the concentra- 
tion parameter. Their triaxial description of the dark halos will be particularly 
useful in predicting a variety of nonsphericity effects, to a reasonably reliable 
degree, including the weak and strong lens statistics, the orbital evolution of 
galactic satellites and triaxiality of galactic halos, and the non-linear clustering 
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of dark matter. Some applications of the triaxial halo model is now in progress 
(Lee & Suto 2002). 

CL3 
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Figure 5. An example of projected particle distribution for a cluster- 
size simulated halo labelled as CL3. The size of each box is Ir^-^ of 
the halo, and particles in the isodensity shells with ps/pcnt = 100, 
2500, and 6.25 x 10^ are plotted on the xy, yz and zx planes (from left 
to right). The bottom panels show the triaxial fits to five isodensity 
surfaces projected on those planes (Jing &; Suto 2002). 

5. Clustering of dark halo on the light-cone 

All cosmological observations are carried out on a light-cone, the null hyper sur- 
face of an observer at z = 0, and not on any constant-time hypersurface. Thus 
clustering amplitude and shape of objects should naturally evolve even within 
the survey volume of a given observational catalogue. Unless restricting the ob- 
jects at a narrow bin of z at the expense of the statistical significance, the proper 
understanding of the data requires a theoretical model to take account of the 
average over the light cone (Matsubara, Suto, & Szapudi 1997; Mataresse et al. 
1997; Moscardini ct al. 1998; Nakamura, Matsubara, & Suto 1998; Yamamoto 
&: Suto 1999; Suto et al. 1999). We take account of all relevant physical effects 
and construct an empirical model for the two-point correlation functions of dark 
halos on the light-cone (Hamana, Colombi, Sz Suto 2001a; Hamana, Yoshida, 
Suto k Evrard 2001b). 
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Figure 6. Projected distribution of dark matter halos with Mhaio > 
Mmin in the hght-cone output of the Hubble volume simulation; Mmin = 
1 X IO^^/j-^Mq, 5 X IO^^/i-^Mq and 2.2 x IO^^/i-^Mq from top to 
bottom. 




Figure 7. Redshift distribution of the cumulative mass function of 
halos identified in the light-cone output of the Hubble volume simula- 
tion (histograms). 
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Figure 6 plots the projected distribution of identified halos from "light- 
cone output" of the Hubble Volume ACDM simulation (Evrard et al. 2002) 
with = 0.04, OcDM = 0.26, ag = 0.9, Qa = 0.7 and h = 0.7. 

Figure 7 displays the redshift distribution of the cumulative mass func- 
tions of the identified halos with M > 2.2 x 10^^/i"^Mo, 5.0 x W^^h'^MQ, and 
1.0 X lO^'^/i~^M0. The total number of the corresponding halos at < z < 3 is 
21,090, 5,554 and 1,543, respectively. For comparison, we also plot the predic- 
tions on the basis of the Press-Schechter mass function {short-dashed lines), the 
modified mass function by Sheth & Tormen (1999; dotted lines), and the fitting 
model by Jenkins et al. (2001; long-dashed lines). The Press-Schechter model 
underpredicts the halo abundance at z > 1, while the Sheth- Tormen model over- 
predicts beyond z ~ 1.5. This tendency is consistent with the previous finding 
of Jenkins et al. (2001) that the Sheth- Tormen model overestimates the number 
of halos when ln{a^^) becomes large. 




1 10 lOM 10 102 

R [h-'Mpc] R [h-'Mpc] 



Figure 8. Two-point correlation functions of halos on the light-cone; 
simulation results (symbols; open circles and filled triangles for real 
and redshift spaces, respectively) and our predictions (solid and dotted 
lines for real and redshift spaces, respectively). The error bars denote 
the standard deviation computed from 200 random re-samplings of the 
bootstrap method. The amplitudes of ^^'^ for Mhaio > 5.0x lO^^/j-^M© 
are increased by an order of magnitude for clarity. 

Figure 8 compares our model predictions with the clustering of simulated 
halos. For the dark matter correlation functions, our model reproduces the sim- 
ulation data almost perfectly at R > 3/i~^Mpc (see also Hamana et al. 2001a). 
This scale corresponds to the mean particle separation of this particular simu- 
lation, and thus the current simulation systematically underestimates the real 
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clustering below this scale especially at z > 0.5. Our model and siniTilation data 
also show quite good agreement for dark halos at scales larger than 5/i~^Mpc. 
Below that scale, they start to deviate slightly in a complicated fashion de- 
pending on the mass of halo and the redshift range. This discrepancy may be 
ascribed to both the numerical limitations of the current simulations and our 
rather simplified model for the halo biasing. Nevertheless the clustering of clus- 
ters on scales below 5/i~^Mpc is difficult to determine observationally anyway, 
and our model predictions differ from the simulation data only by ~ 20 percent 
at most. Therefore we conclude that in practice our empirical model provides a 
successful description of halo clustering on the light-cone. 



6. What's next ? : dark halos vs. galaxy clusters 



Abell clusters 

selected by eyes 
from the photographic plates 
galaxies satisfying 
m,< m < m,+ 2 



T 



SZ clusters 

inverse Compton of CMB 

^Isz~neTeRcl 



Press-Schechter halos 

nonlinear spherical coUapse 



A . = 18 Jt' 

vir 



Halos in N-body simulations 

friend-of-friend method 
linking length 

= 0.2 (mean pailicle separalion) 



1 



X-ray clusters 

thermal bremsstrahlung 



Figure 9. Dark halos - galaxy clusters connection. 



As I have briefly shown in the above sections, understanding of the density 
profiles and clustering of dark halos has been significantly advanced during last 
several years. With those theoretical/empirical successes in mind, the next nat- 
ural question is how to apply them for the description of real galaxy clusters. 
In principle this may be fairly straightforward, but in reality the main obstacle 
is the lack of the proper definition of clusters. This point of view is clearly 
illustrated in Figure 9 which summarizes the conventional definitions of clus- 
ters/halos in various situations. There are a wide range of practical (and quite 
different!) definitions of dark halos and clusters of galaxies. Of course they are 
closely related, but the one-to-one correspondence is unlikely and nothing but 
a working hypothesis. We we need more quantitative justification or modifica- 
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tion of the working hypothesis in order to move on to precision cosmology with 
clusters. 

This problem has not been considered seriously probably because the agree- 
ment between model predictions and available observations seems already satis- 
factory. In fact, since current viable cosmological models are specified by a set 
of many adjustable parameters, the agreement does not necessarily justify the 
underlying assumption. Thus it is dangerous to stop doubting the unjustified 
assumption because of the (apparent) success. 
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